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Abstract 

The proton and neutron density distributions, one- and two-neutron separation 
energies and radii of nuclei for which neutron halos are experimentally observed, 
are calculated using the self-consistent Hartree-Fock-Bogoliubov method with the 
effective interaction of Gogny. Halo factors are evaluated assuming hydrogen-like 
antiproton wave functions. The factors agree well with experimental data. They 
are close to those obtained with Skyrme forces and with the relativistic mean field 
approach. 

1 Introduction 

In the last years, more and more experimental evidence concerning neutron distributions 
in nuclei has become available In particular, it has been found that several neutron 
rich nuclei display a neutron skin or even a neutron 'stratosphere' called the neutron 
halo. Such properties of neutron distributions have been predicted by the asymptotic 
density model 0, the relativistic mean field theory ||, and Hartree-Fock calculations 
with Skyrme forces Q. 

With the increasing amount of experimental data available, halo factors can now 
be used as an additional test of nuclear effective interactions. Namely, microscopic ap- 
proaches should be able to reproduce not only binding energies and charge radii, but 
also the detail of density distributions, especially density tails at large distance from the 
nucleus center. 

*The work is partly supported by the Polish Committee of Scientific Research under Contract No. 
2P03B Oil 12 
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The aim of this work is to perform such a test for the effective interaction proposed 
by Gogny ||J]. While many applications have shown the high quality of this force for 
describing a wide range of nuclear properties [4-7], no systematic study in this direction 
has been made up to now. The mean-field theory employed with the Gogny force is 
applied to nineteen halo nuclei found in |IJ. We calculate neutron and proton density 
distributions at large distance, one- and two-neutron separation energies, and the halo 
factors related to the probability of annihilation of antiproton from atomic-like orbitals. 
These results are compared with experimental data and with Hartree-Fock calculations 
performed with the SLy4 interaction ||. 

In Section 2, we recall the form of the Gogny effective interaction and we describe 
the method by which the results mentioned above are obtained numerically In Section 
3, neutron separation energies and density distributions at large distance from the center 
of the nucleus are analyzed. In Section 4, the prescription employed for evaluating halo 
factors is given and the theoretical values obtained are compared with experimental data. 
Conclusions and plans for further investigations are presented at the end of the paper. 

2 Description of the approach 

The Gogny two-body effective nuclear interaction has the following form H : 
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The first line represents two finite range terms (i = 1,2), with the usual superposition 
of Wigner, Bartlett, Heisenberg and Majorana spin-isospin contributions. r*j is the space 
coordinate of nucleon i, and P a and P T are the exchange operators of spin and isospin 
variables, respectively. The second line of eq. (1) describes a two-body zero-range density- 
dependent interaction. The last line contains a two-body zero-range spin-orbit term and 
the Coulomb potential between protons. Here, V 12 = Vi — V 2 , and 0{ is twice the spin 
operator of nucleon i. 

The set of parameters adopted since 1983, called D1S M, is: 



pLi = 0.7 fm (jl 2 = 1-2 fm 

Wt = -1720.3 MeV W 2 = 103.639 MeV 

B x = 1300 MeV B 2 = -163.483 MeV 

#i = -1813.53 MeV H 2 = 162.812 MeV 

Mi = 1397.60 MeV M 2 = -223.934 MeV 

t 3 = 1390.60 MeV im 3( - 1+ ^ x = 1 

7 = 1/3 W LS = 130 MeV fm 5 
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Hartree-Fock-Bogoliubov (HFB) self-consistent calculations have been carried out with 
this interaction, with the purpose of extracting binding energies and nucleon density dis- 
tributions. Similar results employing the Hartree-Fock (HF) method, i.e. ignoring pairing 
correlations, have also been derived in order to compare them with those obtained with 
the Skyrme HF method. All calculations have been performed in spherical symmetry. The 
method used is described in detail in Ref. ||. The assumption of a spherical mean-field 
may seem inappropriate for those halo nuclei which ground states are deformed. However, 
since halo factors depend essentially on the slope of the density at very large distance from 
the nuclear surface, such a simplification should be sufficient for a first estimate of their 
magnitudes. In the same spirit, the effect on nuclear densities of correlations beyond the 
mean-field approximation, as those coming from oscillations around a fixed shape or from 
large amplitude collective vibrations 0, [TO] has been ignored in the present study. One 



should nevertheless keep in mind that neglecting the increase of neutron and proton radii 
caused by deformation and correlations should lead to a slight underestimation of halo 
factors in most nuclei. 

Our calculations are performed by expanding the HFB one quasi-particle states on 
finite harmonic oscillator (HO) bases. A crucial point in this method is to carefully 
choose the two parameters: number of shells (Nmax) and harmonic oscillator frequency 
hu), on which these bases depend. This is especially needed in the present study where the 
behaviour of single-particle wave-functions at large distance from the center of the nucleus 
has to be accurately determined. The method we have employed consists of choosing, for 
each value of Nmax, the hu value that minimizes the HFB binding energy. Then Nmax 
is increased until convergence of the HFB binding energy is obtained. 

Application of this method to 58 Ni is illustrated in Figs. 1, 2 and 3. Fig.l displays the 
variations of the HFB energy (B) with hu for Nmax=16. The energies of the minima of the 
curves B{hu) are plotted in Fig.2 as functions of Nmax- One observes that convergence 
of the HFB binding energy is obtained for Nmax=16 in this nucleus. As shown in Fig.3, a 
well defined energy minimum is then obtained in the (hu, N) plane. This procedure has 
been employed for the nineteen halo-nuclei considered in the present study. 

In order to illustrate the influence of the HO basis parameters on densities, we have 
plotted in Fig. 4 the logarithm of the neutron density distribution log p n of 58 Ni versus 
the radial distance r for a few values of hu. One can see that p n is sensitive to the chosen 
hu only in the peripheral region r > 8 fm, a region where p n is smaller than 10 -10 fm~ 3 . 



3 Results 

The parameters of the HO bases employed for all the halo-nuclei studied in the present 
work are gathered in Table 1. The HO frequency hu corresponding to the minimal HFB 
energy and the number of shells Nmax above which the energy does not change any more 
are shown for each nucleus. The HFB energy and the one- and two-neutron separation 
energies S n and S2 n are also listed. The separation energies have been found by subtraction 
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of the HFB energies of neighbouring isotopes 



S n (Z, N) = B(Z, N) - B(Z, N-l), (2) 

S 2n (Z, N) = B(Z, N) - B(Z, N - 2) , (3) 

Odd isotope energies have been calculated using the blocking version of the HFB theory, 
as described in Ref . [§J . In each nucleus, the blocked quasi-particle state has been chosen 
as the one having the experimentally known spin nearest to the neutron Fermi surface. 



Quadrupole deformations fa taken from Ref. [ I L| are given in Table 1 in order to indicate 
which nuclei are deformed in their ground states. Let us recall that all nuclei are considered 
as spherical in our calculation. The last two columns of Table 1 display calculated and 
experimental halo factors, the meaning of which will be explained in Section 4. 

The one- and two-neutron separation energies of these nuclei are compared with 



the experimental values [12] in Figs. 5 and 6. In these figures, the differences between 
experimental and theoretical separation energies are plotted as functions of the proton 
number Z. 

Calculated neutron separation energies agree with experiment within ±1 MeV in spher- 
ical nuclei, except for 48 Ca, 96 Ru, 130 Te and 144 Sm. In these nuclei, the neutron Fermi 
level is located in the vicinity of a major shell gap, where the single-particle level density 
is underestimated in the present mean-field approach. A proper description of separation 
energies in nuclei near closed shells would require to take into account correlations beyond 
the mean-field - as RPA ground state correlations - and the induced spectroscopic factors 
and single-particle level displacements. Let us also note that, in the calculation of S n s, 
odd nuclei have been assumed to be spherical, which leads to a systematic overestimation 
of theoretical values. In the case of deformed nuclei, one observes that S n s usually are 
overestimated, while S^s often are underestimated with the present spherical approach. 

Calculations of proton and neutron density distributions have been performed for 
the nineteen halo nuclei using the HF and HFB procedures with the Gogny interaction. 
For the purpose of comparison, similar results have been derived also with the SLy4 
parametrization || of the Skyrme force ]13) using the HF method. 

The proton and neutron density distributions p n ( p ) for 58 Ni are presented in Fig. 7. 
The Skyrme interaction gives proton and neutron densities lower by ~ 10% in the nucleus 
interior and, consequently slightly larger proton and neutron radii than the Gogny force. 
However, the tails of the densities obtained within the two approaches do not differ sig- 
nificantly from each other. Let us recall here that the parameters of the Gogny force have 
been determined in order to allow for the inclusion of ground state correlations in the 
description of one-body observables. As a consequence, nuclear radii are expected to be 
accurately reproduced by the Gogny force only when correlations beyond the mean-field 
are included. 

The difference between the proton and the neutron distributions at large distance is 
important for the investigation of halo effects. In order to illustrate the detailed structure 
of these distributions, several functions of the densities p n (p) which enter halo factors (see 
Section 4) have been plotted. A complete set of such results is displayed in Fig. 8 for 
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160 Gd as an example. Similar plots for the remaining eighteen nuclei are shown in the 
next eight figures. 

The upper part of Fig. 8 represents the density distributions for protons (solid lines) 
and neutrons (dashed lines) obtained with the Gogny force without pairing correlations 
(GoHF), with pairing correlations included (GoHFB) and with the Skyrme force SLy4 
(SkHF) without pairing correlations. The densities are plotted up to r = 12 fm, a ra- 
dial distance above which they are lower than 10~ 10 fm~ 3 . The middle row shows the 
logarithms of these densities, and the leftmost diagram of the lowest row the difference 
logpn - log Pp. 

One can see that the densities computed with and without pairing correlations are 
very close to each other, although the contribution of pairing to the binding energy is of 
the order of 10 MeV. 

More significant differences can be noticed between the results obtained with the SLy4 
and the Gogny forces (both without pairing), especially for large r and near r = 0. 

It is interesting to analyze how the different single-particle states contribute to the 
total density p. Denoting by p u the contribution of single-particle state v : 

P = E^ (4) 

V 

the relative contributions p v jp of the occupied single-particle orbits are plotted in the 
rightmost diagram of the lowest row in Fig. 8. These results have been obtained from 
HF calculations with the Gogny force. One observes that all single particle states except 
one, contribute more or less the same amount in the whole region r = — 12 fm. At 
large distance and for r ~ 0, the f?/2 one-neutron state strongly dominates all the other 
ones. This clearly indicates that the halo is a single particle effect in this nucleus, which 
confirms earlier results obtained with Skyrme forces [Q] and with the relativistic mean 
field theory §. 

The middle plot of the lowest row in Fig. 8 shows the function (p n — p p ) ■ r 2 / (N — Z) 
calculated within the three approaches GoHFB (solid line), GoHF (dashed line), and 
SkHF (dotted line). This function is strongly correlated with the magnitude of halo 
factors. 

Similar results for the other eighteen halo nuclei listed in Table 1 are shown in the 
next eight figures: 9a, 9b, 10a, 10b, 11a, lib, 12a, and 12b. Figures "a" correspond to 
the lighter Ca-Sn nuclei and figures "b" to the heavier Sn-U. Every multiplot presents the 
same quantity versus r for nine different nuclei. 

In Fig. 9, the differences logp n — logp p between the logarithms of the neutron and 
proton densities are shown. The solid lines have been obtained with the Gogny force 
(GoHF), and the dashed ones with the SLy4 (SkHF) interaction. The difference between 
the proton and neutron densities grows with r and is larger for the Skyrme force in lighter 
nuclei (9a) while, for the heavier ones (9b), the Gogny force gives similar or larger density 
differences than SLy4. 

Fig. 10 shows the contributions p v of neutron (solid lines) and proton (dashed lines) 
occupied single particle orbitals v to the whole density p. One can see that in most cases 
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only one neutron state determines the magnitude of the density tail. For lighter nuclei 
(10a), only 116 Cd and 100 Mo have more than one neutron state contributing to the tail of 
the density. In heavier nuclei, Fig. 10b, the Sn and Te isotopes show a similar behaviour, 
with two neutron orbitals contributing to the density tail. A very interesting situation 
is found in 144 Sm where the contribution that dominates corresponds to a proton state. 
Therefore a proton halo can be expected in this nucleus. This is in line with experimental 
data |T[. Unfortunately this is not the case for 106 Cd, which also has a proton-rich nuclear 
stratosphere. For every plot in Fig. 10, the quantum numbers Ij of the orbital for 
which the contribution p v is maximum at large distance is indicated. The single particle 
character of the nuclear periphery has also been found in HF calculations with the Skyrme 
force Skill [l[ and within the relativistic mean field theory ||. 

The two Fig. lla,b show the function (p n — p p )r 2 /(N — Z) versus r, for the eighteen 
halo nuclei. This quantity directly enters the integral for the halo factor, Eq. (6), and 
determines the neutron or proton halo in the nuclear periphery. HF results are shown 
both for the Gogny (solid line) and the SLy4 (dashed line) forces. It is apparent on 
these curves that the Skyrme force yields a slightly larger neutron halo than the Gogny 
interaction. 

Figs. 12a and 12b illustrate the nuclear skin effect. The r 2 -weighted difference between 
the average neutron (p n /N) and proton (p p /Z) single particle densities calculated with 
the Gogny force is plotted for each nucleus as a functions of r. One can see that this 
quantity strongly oscillates in the surface region, with variations that significantly differ 
in amplitude and shape, depending on the nucleus considered. 



4 Halo factor 

In experiments probing the nuclear periphery using the formation of antiproton atoms, 
antiprotons are catched on hydrogen-like atomic orbitals and then annihilated on the 
outer tail of the proton or neutron density. The halo factor is defined as |j] : 

s 

where the summation goes over all the antiprotonic states with meaningful annihilation 
widths r^j on proton {p) or neutron {n) and s denotes the atomic state of the antiproton. 

The widths are calculated by integrating the density of one kind of nucleon with the 
square of the antiproton wave-function ip s (r) : 

K( P ) = J Pn( P M ip) (r)\ 2 P(r)r 2 dr. (6) 

tp s (r) is taken as the solution of the Schrodinger equation for a hydrogen-like antiprotonic 
atom. The nucleus is assumed to be spherical. Let us point out that the nuclear radius 
is much smaller than typical antiproton orbital radii. 
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The factor P(r) describes the probability that the annihilation products of the an- 
tiproton will not be absorbed by the nucleus. This factor, which should depend on the 
atomic state of the antiproton, will be assumed here to be independent of it. P(r) takes 
into account the deep hole creation probability Pdh(^) and the pion escaping probability 
-P7r,esc(f )• For the ground state nuclear periphery, only cold annihilation is important. Hot 
annihilation has to be eliminated from the widths T, which leads to the P(r) = constant 
approximation. 

In Fig. 13 the halo factors obtained with the Gogny and Skyrme forces are com- 
pared with the experimental data taken from Ref. The variations of experimental 
halo factors from nucleus to nucleus are satisfactorily reproduced by the two theoretical 
calculations. In most cases, the Gogny force gives smaller halo factors than the Skyrme 
interaction, which is consistent with the fact that the Gogny D1S parameterization under- 
estimates radii when correlations beyond the mean field are neglected. This is especially 
true in the five heavier nuclei which all are deformed. The halo factors derived with 
SLy4 appear in good quantitative agreement with experimental data. In particular, the 
unusually large halo factor in 176 Yb is well reproduced. Let us note that both effective 
interactions significantly underestimate the halo factors of the two Te isotopes, which 
indicates that the halo structure of these two spherical nuclei is not properly described 
by the present microscopic calculations. One reason for these discrepancies may be that 
correlations beyond the mean-field approximation are necessary for a correct description 
of neutron densities at large distance in these nuclei. 

5 Conclusions 

Hartree-Fock-Bogoliubov calculation in spherical symmetry have been performed with 
the Gogny force for nineteen halo nuclei. The parameters of the HO bases employed - 
maximum number of shells Nmax and frequency hu - have been carefully chosen in order 
to ensure the convergence of the total binding energy, and to describe the large distance 
behaviour of nucleon density distributions. 

The following conclusions can be drawn from the present investigation: 

1. One- and two- neutron separation energies are reproduced with the Gogny force 
within ± 1 MeV in most of the spherical nuclei studied. In deformed nuclei, the 
spherical symmetry assumption leads to an underestimation of S-2 n s and an overes- 
timation of S n s that can reach 3 MeV. 

2. The Gogny force D1S gives an overall satisfactory account of the nuclear periphery 
which could probably be improved by taking into account deformation effects and, 
to a lesser extent, correlations beyond the mean field. Inclusion of pairing correla- 
tions in the self-consistent calculation induces only weak changes in the proton and 
neutron density distributions, and has a negligible effect on halo factors. 

3. The variations of halo factors from nucleus to nucleus obtained with the Gogny force 
and the SLy4 interaction are in good agreement with experimental data. From the 
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present spherical mean-field calculations, SLy4 gives a larger neutron halo effect 
than the Gogny force, in better agreement with experiments. 

As almost half of experimentally known halo-nuclei are well deformed, a test of the 
influence of including nuclear deformation in the self-consistent calculations is clearly 
needed. In the case of the Gogny force, a similar test concerning the effect of the ground 
state correlations associated with oscillations of the mean-field - large amplitude collective 
motion in soft nuclei, RPA ground state correlations in rigid nuclei - and of possible shape 
coexistence phenomena should also be performed. This is left for future work. 
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Table captions: 



1. The results of the Hartree-Fock-Bogoliubov calculation with the Gogny D1S force 
for 19 halo nuclei [l]. The harmonic oscillator basis parameters: the shell number 
N Mj 4x and the oscillator frequency hu, the binding energies B, one S n and two S 2n 
neutron separation energies, the equilibrium deformation j3 2 and the theoretical / 
and experimental halo factors are listed in the table. 

Figure captions: 



1. Hartree-Fock energy (crosses) obtained with the Gogny force in 58 Ni versus the fre- 
quency (hu) of the harmonic oscillator basis. The theoretical points are interpolated 
by 2 nd order polynomial in Tiuo. 

2. Hartree-Fock energy (crosses) obtained with the Gogny force in 58 Ni as a function 
of the number (Nmax) of harmonic oscillator shells included in the basis. 

3. Contours of the Hartree-Fock energy (crosses) obtained with the Gogny force in 58 Ni 
in the (N M ax, ^) plane. N M ax is the number of shells included in the harmonic 
oscillator basis and hu is the harmonic oscillator frequency. 

4. Logarithms of the neutron density distribution p n as functions of the radial dis- 
tance r in 58 Ni. The different curves correspond to different values of the harmonic 
oscillator basis frequency fiui. 



5. Difference between theoretical (Gogny) and experimental [|1^] one-neutron separa- 
tion energies S n for experimentally known halo-nuclei. 



6. Same as Fig. 5 for two-neutron separation energies 5": 



2n- 



7. Proton (dashed and dotted lines) and neutron (solid, thin dashed lines) density 
distributions obtained from Hartree-Fock calculations with the Gogny force (GoHF) 
and with the Skyrme SLy4 force (SkHF), versus the radial distance r. 

8. Density distribution results obtained in 160 Gd. First row: Density distributions for 
protons p p (solid lines) and neutrons p n (dashes lines) obtained from Hartree-Fock 
calculations with the Gogny force D1S (first column), from Hartree-Fock-Bogoliubov 
calculations with the Gogny force (second column), and from Hartree-Fock calcula- 
tions with the Skyrme force SLy4 (third column). Second row: Logarithms logp n ( p ) 
of the densities shown in the first row. Third row, first column: the difference 
logp n — logpp between the logarithms of the neutron and proton densities; sec- 
ond column: the function (p n — p p )r 2 /(N — Z) obtained by applying the Hartree- 
Fock method with the Gogny force D1S (GoHF, dashed line), the Hartree-Fock- 
Bogoliubov method with the Gogny force (GoHFB, solid line) and the Hartree-Fock 
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method with the Skyrme force SLy4 (SkHF, dotted line); third column: single- 
particle contributions p v to the density p for protons (solid lines) and neutrons 
(dashed lines). The Ij quantum numbers of the orbit with largest p v at r = 12 fm 
are indicated in each plot. 

9a. The differences logp n — logp p obtained from Hartree-Fock calculations with the 
Gogny force D1S (GoHF - solid lines), and with the Skyrme force SLy4 (SkHF - 
dashed lines) for the nine lighter halo-nuclei 48 Ca to 112 Sn. 

9b. Same as Fig. 9a for the nine heavier halo-nuclei 124 Sn to 238 U. 

10a. Relative contributions p v j p of the single particle proton (dashed lines) and neutron 
(solid lines) states v to the total density p as functions of the radial distance r for 
the nine lighter halo-nuclei 48 Ca to 112 Sn. 

10b. Same as Fig. 10a for the nine heavier halo-nuclei 124 Sn to 238 U. 

11a. The function (p n — p p )r 2 / (N — Z) obtained from Hartree-Fock calculations with the 
Gogny force D1S (solid lines), and with the Skyrme force SLy4 (dashed lines) versus 
the radial distance r for the nine lighter halo-nuclei 48 Ca to 112 Sn. 

lib. Same as Fig. 11a for the nine heavier halo-nuclei 124 Sn to 238 U. 

12a. The function (p n /N — Pp/Z) r 2 obtained from Hartree-Fock-Bogoliubov calculations 
with the Gogny force for the nine lighter halo-nuclei 48 Ca to 112 Sn. 

12b. Same as Fig 12a for the nine heavier halo-nuclei 124 Sn to 238 U. 

13. The halo factors of experimentally known halo-nuclei obtained from Hartree-Fock- 
Bogoliubov calculations with the Gogny force D1S (stars joined by the thin dashed 
lines) and from Hartree-Fock calculations with the Skyrme forces SLy4 (crosses 
joined by dashed lines), compared with experimental data (plusses wit errorbars 
joined by solid lines) from Ref. [1]. 
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